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We resolve an apparent contradiction between numeric and analytic results for one-dimensional 
disordered systems with power-law spectral correlations. The conflict arises when considering rig- 
orous results that constrain the set of correlation functions yielding metallic states to those with 
non-zero values in the thermodynamic limit. By analyzing the scaling law for a model correlated 
disorder that produces a mobility edge, we show that no contradiction exists as the correlation 
function exhibits strong anticorrelations in the thermodynamic limit. Moreover, the associated scal- 
ing function reveals a size-dependent correlation with a smoothening of disorder ampHtudes as the 
system size increases. 



I. INTRODUCTION 

The necessary and sufficient conditions for the occur- 
rence of a disordered Anderson transition have received 
a fair amount of attention in the past half century.'^ltSl 
Early theoretical work with uncorrelated and short-range 
correlated disorder potentials identified /3 = ^}^^ as 
the relevant scaling function. If (3 depends solely on 
the dimensionless conductance G, all the eigenstates of 
a one-dimensional system will be localized.!^ This as- 
sumption is known as the single parameter scaling (SPS) 
hypothesis.l^ It was realized early on that a rigorous 
proof of the SPS hypothesis involves an analysis of the 
full probability distribution of InQ, or equivalently, its 
cumulants.^*^ For one-dimensional systems (with uncor- 
related disorder potentials), an analysis of these quan- 
tities helped determine a relation between the average 
dimensionless conductance and its variance, V ar (In G) — 
-2(lng), valid when the SPS hypothesis holds .li^liU This 
relationship fails, howe ver^ fo r states at the band edge,'^ 
and at the band centei'^^lH! which in turn signifies the 
failure of the SPS hypothesiiSl at these energies. These 
violations allow the possibility of extended states but do 
not guarantee the presence or absence of a delocaliza- 
tion transition (mobility edge) . More complications arise 
in the analysis of the SPS hypothesis for systems with 
correlated disorder where violations of SPS appear to be 
correlated with a crossover between two different scaling 
regimes for the localization length.'^^Hi^ 

Due to these limitations, different authors have en- 
deavored to determine the necessary conditions for the 
existence of an Anderson transition from a more rigorous 
perspective. A seminal work was carrie d ou t by Kotani 
in his theory of random ergodic operators.'^^In it, he was 
able to show that for the existence of an Anderson transi- 
tion, Gaussian disorder potentials with correlations must 
be deterministic (the whole disordered system can be de- 
scribed from the behavior of the potential within a small 
region) .I^MSl Kotani's theorems also imply that a metal- 
lic band can not exist if T{x), the correlation function 



for the disordered potential, goes to zero as the distance 
x — > oo as a power law or faster.fSi The theorem thus pro- 
vides the mathematical justification for the absence of 
a transition in models with uncorrelat ed and scale-free 
disorder distributions in one-dimension.^^ l^^ l ^^ l More re- 
cently, a sufficient condition was also proposecP for con- 
tinuum disordered models: to produce a mobility edge, 
the disorder potential must be V{x) £ , with /? > 1/2 
and representing the class of continuous functions 
that are /3-differentiable. Along with these theoretical 
advances, new experimental studies on ultra-cold atom 
systems have called into attention the conditions for the 
existence of an Anderson localization transition in the 
presence of correlations.'^^^ 

Several numerical works on ID systems with corre- 
lated potentials have also addressed the conditions for 
the existence of extended states. A broad classification 
distinguishes three groups of models: discrete,-"'^ quasi- 
periodic,'22l and long-range spectral correlations .'^^HUlj^e. 
markably, the latter, characterized by the use of a dis- 
order spectral density S{k) ~ 1/fc", (with a a measure 
of the range of correlations) has been shown to produce 
a localization/delocalization transition. Using a tight- 
binding model in one dimension, a mobility edge is pre- 
dicted to appear for disorder strengths W < At and 
a > 2 (with t the hopping energy). Despite consid- 
erable efforts however, the exact phase diagram for this 
model for W > At remains controversial with contra- 
dictory evidenc^^mH regarding the region of parameter 
space where the transition should occur. Furthermore, 
numerical studies have remarked the smoothening of the 
disorder amplitude generated in this model as the ther- 
modynamic limit is reached, suggesting the transition to 
be order/disorder instead of a true Anderson transition.!^ 

In addition to these issues, the authors in Ref. [7] 
pointed out an apparent disagreement between Kotani's 
theorem and the localization/delocalization transition 
found in the model with ~ l//c" correlations. The argu- 
ment presented makes use of the analytic Fourier trans- 
form of S{k) to find the corresponding correlation fimc- 
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tion in real space, with the introduction of an appropriate 
cutoff in order to deal with the non-analyticity at fc = 0. 
After straightforward algebra, it is found that for values 
of a < 1, the correlation function decays as l/x^~°' . For 
a > 1, however, the Fourier transform renders exponen- 
tially decaying correlation^ as 



r(a;) « e-l'^l' for x>a 



(1) 



where c and b are parameters and a is the short-range 
cutoff. Note that the correlation function goes to zero 
in the thermodynamic limit in both regimes, i.e., r(a;) 
appears to violate the necessary conditions for the exis- 
tence of a mobility edge.^^ The authors thus raised the 
question of the validity of the Anderson transition found 
numerically in Ref. i23i 

In order to resolve these inconsistancies and provide 
some insight into the nature of the transition produced 
by these correlations, we have carried out a study of the 
properties of the discrete real-space correlation function 
as presented in Ref. [531 In this work we present a detailed 
analytic calculation of the scaling expression for the dis- 
crete form of the correlation function, and support the 
resulting expressions by the numerical evaluation of the 
correlation function as defined in the original work. Our 
findings settle the apparent contradiction between the 
necessary conditions as imposed by Kotani's results and 
those obtained in numerical studies of disordered models 
with iS'(fc) ~ spectral functions. We show that the 

model contains anticorrelations between the energies of 
the two most distant sites in the thermodynamic limit. 
We also confirm that the correlation function obtained 
in real space is system-size dependent and exhibits an 
unphysical smoothening of disorder amplitudes as the 
thermodynamic limit is reached. We argue that these 
features, some already identified in previous workSjf^Sl ap- 
pear to play a significant role in the origin of extended 
states. 



colored noise models with a power spectrum P{ujk) = 
i-o^"- Here ujk = kAuj is a discrete series of frequencies 
labelled by the integer k, Auj = 2tt/T, T being the total 
time interval in which the random function is calculated 
as a discrete time series.^* Within this approach, a 
random variable X{tj) {tj is time and j an integer), is 
defined by the discrete Fourier transform 



N/2 

^ cos{uJktj) ■ 



bksm{ujktj) 



(3) 



where N — T/At, At — (tj+i—tj), and the Fourier coeffi- 
cients are defined by the power spectral function as: = 
^/PiuJk)Auj cos{(j)k) and = -^P(cjfc)Aw sm{(j)k)- No- 
tice that (f)k is a random independent variable chosen in 
the interval [0,27r]. Replacing the expressions for the 
Fourier coefficients above yields: 



N/2 



^ ^\/i^k''^^^'^^i^ktj + (kk)- 



(4) 



k=l 



A map to a tight-binding model in real space is straight- 
forward using the transformations: Xj — > tj and recipro- 
cal variable — >■ lo^ 



N/2 

E 

k = l 



Kf. "AKk COs{KkXj + 4>k)- 



(5) 



In this case, N , the numbers of total time-steps is re- 
placed by the total number of lattice sites N = L/a with 
L the system size, a the lattice constant, k^. = 27rfc/L, 
An = 27r/L, and Xj = ja. Replacing these values into 
the above equation we arrive at 
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k=l 



L 



(1-q)/2 



-q/2, 



II. SCALING FORM OF THE CORRELATION 
FUNCTION 

Let us review the procedure to numerically generate 
random on-site correlated variables for tight-binding An- 
derson Hamiltonians. To set up notation, we intro- 
duce the one-dimensional tight-binding model Hamilto- 
nian with N sites as: 



N 

c\cj +h.c. + ^e,c\ci 



(2) 



<i,j> 



where t is the hopping integral, c\ (cj is the creation 
(destruction) operator for one fermion at site i and 
is the on-site disorder energy. < • • • > stands for sum 
over nearest neighbors sites only. In order to generate 
the discrete random variables e^, we follow a standard 
procedure, as presented in Ref. |28| which is based on 



where we have used to emphasize that this is the on- 
site disorder energy (random variable) for the Anderson 
tight-binding model. The disorder average of this quan- 
tity is given by 



N/2 

E 



27r 

TV 



(l-a)/2 



-q/2/ 



(cos(^+0,)) (7) 



N/2 

E| 

k=l 



27r 



(l-a)/2 



((cOS(/)fe) 



COS — (sm (t>k) sm j . 



(8) 



The expectation values {cos <j>k) = (sin^j.) — resulting 



in (ej) = 0. 
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We can calculate the covariance through similar means 

JV/2 W/2 



EE 

q=l k=l 



|27r 



-a/2. 



(cos (^^^ + 0fc) COS + 



(9) 



For k q, the expectation value can be factorized as 
(j)k and (j)q are independent variables and the only term 
contributing to the expression is fc = q: 

N/2 



|27r 
I TV 

k=l 



(10) 



X (cos ^ + cos (^^^ + 
which can be simplified to yield a covariance of 

N/2 



fc=i 



N 



2'Knk 
N ■ 



(11) 



This even function depends solely on the difference be- 
tween positions n ~ \m — j\. The normalized correlation, 
or auto-correlation, function is defined as (for zero mean) 



r(m,j) 



which when combined with Eq. ( 11 ) gives 

Efi?fc""C0s2-"'= 



r(Q;, n) 



N 



(12) 



(13) 



From this point on we impose periodic boundary con- 
ditions with the correlation function defined in n e 
[0, A^/2]. This expression has a natural limit of — ^ oo 
when written in terms of the variable 7 = 2n/N defined 
in [0, 1]. The introduction of this variable allows for the 
explicit removal of the system size dependence in the ar- 
gument of both sums and emphasizes the scaling form 
for the correlation functionPS! 

The thermodynamic limit is obtained by taking the 
upper limit of the sum (N) to infinity. 



r(a,7) 



(14) 



In order to study this expression, we consider two cases: 
a > 1 and a < 1. When a > 1 both the numerator 



and denominator are well defined and Eq. ( 14 ) can be 



rewritten in terms of imaginary exponential functions: 



Y^oo l.-a( iirf\k ,-a( -iTTf\k 

r(a,7)^ ^^-^' 2Erifc- ^''^ 

The terms in the numerator are the polylogarithm func- 
tions Lia{z) while the denominator corresponds to the 




FIG. 1. Plot of the correlation function as a function of the di- 
mensionless distance 7. The red solid line corresponds to the 
critical value of a = 2 which is given by Eq. | |19[ ). The other 
curves correspond to various values of a in Eq. (16 1. (In- 



set) Correlation function between the two most distant points 
(7 = 1). The origin corresponds to the critical value a = 2 
and yields a correlation of —1/2. For larger values of a, the 
correlation becomes more negative. 



Riemann-Zeta sum C(«)- The final expression for a > 1 
is then 



r(a,7) ^ 



£»a(e^^^) +£»a(e~^^^) 
2C(a) 



(16) 



For general values of a > 1, Eq. (16) cannot be sim 



plified further, however, it takes a very simple expression 
at a = 2. It is important to remark that this value corre- 
sponds to the critical value beyond which a band of ex- 
tended states appears.!^ By introducing the second order 
Bernoulli polynomial, ^2(7/2), and using the identity 



(^27^)^ 



^2(7/2) 



we obtain 



(17) 



(18) 



where C(2) = 7rV6 and ^2(7/2) = (7/2)^ - 7/2 + 1/6. 
Finally, the correlation function reduces to: 



r(a-2,7) - 72-37+1. 



(19) 



This expression is plotted in Fig. [T] together with re- 
sults from Eq. ( 16 ) for a few values of a above and below 
the critical value a = 2. Three features are distinguished 
in these results: 

1. ) The function is linear near 7 ^ for the critical 
value of a = 2 while the behaviour is convex for a > 2 
and concave for a < 2. 

2. ) The correlation function goes negative for a > 1. 

3. ) The correlation function converges to a non-zero 
value at the thermodynamic limit of 7 = 1 and yields a 
value of —1/2 at the critical value, a — 2. 
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FIG. 2. Typical disorder realization for N = 1000, a = 2 
(solid, blue), a = 4 (dotted, red), and a = 10 (dashed, black). 
Energy is normalized as (e^) = 1. Plots shows tendency for 
strong correlations to induce a sinusoidal structure in the dis- 
order configuration. 



Notice that features 1.) and 2.) occur for smaller val- 
ues of 7 and as such can be considered as short-range 
effects. Feature 3.) however, corresponds to a non-zero 
value of the correlation function at infinite range. 

Furthermore, the value of the correlation function be- 
tween the two most distant points, 7 = 1, is 



r(a,7 = l) 



C(«) 



(20) 



where Lio,{ — l) ~ —r\(oi) with 77(a) the Dirichlet-eta func- 
tion. We can relate r\{pi) to the Riemann-zeta function 
through the relationship r\{a) = (1 — 2^~")C(q;) and ob- 
tain 



r(a,7 = 1) = 2^- 



1 



(21) 



which is plotted in the inset of Fig. [T] 

The most remarkable characteristic of these expres- 
sions is the existence of a very strong negative correlation 
between infinitely separated sites. These negative corre- 
lations indicate that any two sites separated by distances 
of the order of the system size are statistically more likely 
to have energies with opposite signs. As a result, an 
overarching sinusoidal structure begins to develop in the 
disorder configuration that allows for the appearance of 
extended low energy states when a > 2. We assign this 
structure to be one of the causes for the emergence of 
extended states as illustrated in Fig. [2] which highlight 
the smoothening of potential amplitudes with increased 
system size, already proposed in previous works:^ We 
should note also that a second important element ap- 
pears to be the value of the correlation strength at infi- 
nite distances: our results indicate that it should exceed 
the value of r(7 = 1) = —1/2 to produce extended states. 

We turn our focus to the case when a < 1. We return 
to Eq. ( |14[ ) and examine the convergence of the numer- 
ator and denominator. The denominator is easily seen 



to be divergent as the series converges more slowly than 
the harmonic series for any values of a < 1. The con- 
vergence of the numerator can be seen by applying the 
Dirichlet convergence test for all a > 0. Note that the 
series for a = is not convergent but bounded. Thus, for 
all values of 7 > the correlation function is since the 
numerator is bounded (or convergent) and the denomina- 
tor diverges. At 7 = 0, the numerator and denominator 
are equal and the correlation function gives the value 1. 
The correlation function for a < 1 can be summarized 
by: 



r(a,7) 



37,0 



(22) 



where b-^^ is the Kronecker-delta function. Thus for val- 
ues of a < 1, the correlation function corresponds to an 
effective short-range correlated function that precludes 
the existence of extended states. 

Lastly, an additional assumption necessary for the ap- 
plication of Kotani's results is the fact that the random 
variables generated with correlations given by Eq. ( 13 ) 



must be Gaussian distributed. A straightforward test for 
this condition consists of calculating the fourth cumulant 
of the distribution to check if its value is zero. By follow- 
ing a similar procedure as the one outlined above, one 
can show that the fourth cumulant of the distribution is 
7^4 = |C(a)^ - |C(2a) when a > 1. Note that for a < 1 
the correlations correspond to the regime where no tran- 
sition takes place. Because it has a non-zero value, the 
distribution generated is effectively non-Gaussian thus 
ruling out the applicability of Kotani's theorems. 

Finally, and in addition to the analysis of the scaling 
form, it is important to mention that our study of the 
/3-function and its variance for this model reveals strong 
violations of the SPS relation Vari^Y^Q) = — 2(lnC/) for 
all energies when a > 1. This result further highlights 
the peculiar nature of the transition found for the set of 
values a > 2. 



III. DISCUSSION 

An identification of the salient features of potential 
models with ~ 1/k" correlations as proposed in Ref. US 
shows that the disordered potential thus generated does 
not contradict the rigorous conditions put forward in 
Kotani's work. The scaling form of the correlation func- 
tion does not vanish in the thermodynamic limit (a nec- 
essary condition for the absence of metallic states), but 
reaches negative values. Moreover, the random energies 
are not Gaussian variables as required by the theorem. 
As a consequence, Kotani's theorem can not be used to 
rule out the existence of a mobility edge for this type of 
model. Notice that a key result of our analytic derivation 
in the limit — > 00 is the scaling form for the correlation 
function. This scaling emphasizes the fact that the cor- 
relation between any two given points is size-dependent, 
producing a smoothening of disorder amplitudes as the 
system size increases. This feature, that has been already 



5 



pointed out in previous workJlSEni favors the rise of an 
ordered potential in the thermodynamic hmit, with its 
corresponding extended-Uke states. 

More insight into the effect of long-range negative cor- 
relations can be gained by anal yzi ng the correlation func- 
tion of the random dimer model^with site energies ei, £2- 
As it is well known, a discrete level with extended states 
appears at ei = t and 62 = —t. In this case, a numeri- 
cal evaluation of the correlation between any two points 
yields ^ 0.3. Once again, this value is different from zero 
and consistent with Kotani's theory. 

Finally, these results suggest that two key ingredients 
are necessary for the existence of a band of extended 
states: the crossover to negative correlations in the ther- 
modynamic limit, and a minimum correlation strength 



between the two most distant points. When considered 
together, however, they indicate the transition to be an 
order/disorder one instead of a classical Anderson tran- 
sition. 
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